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Abstract. We introduce the concept of viscosity (both shear and bulk) in the context of hadron 
physics and in particular the meson gas, highlighting the current theoretical efforts to connect 
O^l possible measurements of the viscosities to underlying physics such as a phase transition or the 

trace anomaly. 
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CONCEPT OF VISCOSITY 

The reader might be familiar with the Navier-Stokes equation, easily recognizable as a 
^! Newton's force equation for an element of viscous fluid 

9 ( Tt + (v ' v)v ) = ~ VP+ vAv + + 3 ) v(v ' v) • (1) 

! The constant 77 is the shear viscosity. It controls the flow of the momentum component p y 

perpendicular to the direction of flow v x as per j = — f] 7^ ■ This is a first order transport 
equation, that assumes that the higher derivatives of the velocity are of negligible 
importance. If the fluid is compressible, then a second transport coefficient, the bulk 
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viscosity £, controls the relaxation of longitudinal momentum gradients. 

For a dilute gas, momentum transport is diffusive (a particle flows along carrying 
its momentum) and is hampered by the scattering cross section that interrupts the free 
streaming, for a non-relativistic hard-sphere gas 77 = \nmvX in terms of the mean free 
path X = -L. Thus, the larger the interaction, the smaller the viscosity; in such a diffusive 
system the "perfect fluid" and the "perfect gas" are opposite limits, characterized by 
strong the first and weak interactions the second, by low and high viscosity respectively 
(provided that the hydrodynamic description still applies). 

Maybe it is interesting to note here that the dissipation of mechanical energy into heat 
is directly proportional to these viscosities 
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(a temperature gradient adds a term proportional to the thermal conductivity K, and so 
on for other transport coefficients). 
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To generalize the definition of viscosity to a relativistic system [1] it is more conve- 
nient to look at the momentum-stress tensor than the force Eq. (1). For a viscous fluid, 
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n ; - fc = pVjV k + P8 ik - T]{y i)k + v Ki - -8 ik v l j) - C$ikV j (2) 

that can be carried over to special relativity with the substitution of the density p by the 
enthalpy density w because of the mass-energy equivalence, and the velocity v by the 
four- velocity uq , 

Tap = -Priap + wu a up + T aj6 (3) 

with T a p the dissipative part of the momentum-stress tensor. 

This relativistic generalization suffers from a frame-definition ambiguity (choice of 
four- velocity up) that affects how dissipation is apportioned between thermal conduc- 
tivity and bulk viscosity, for example. We follow Landau's convention that the proper 
frame is that in which the momentum flow vanishes, and the energy satisfies ideal hydro- 
dynamics (no dissipation). This amounts to choosing T a $u$ = and lifts the ambiguity. 
Then the thermal conductivity does not appear in the dissipative part of the stress-energy 
tensor (it is confined to the particle-number dissipative flow) 

= -ri(u a p+up^+upu 7 u ajr + UaU 7 up r )- U-\*\ (lap + u a up)u J \ y . (4) 

The fluid equations of motion are then the conservation law for the stress-energy 
tensor T a ^ y obtained by combining Eqs. (3) and (4). 



THE KSS NUMBER 



In fluid mechanics it is common to construct dimensionless ratios that allow similarity 
analysis between different systems. Perhaps the best known of these is the Reynolds 
number R = that quotients the density, characteristic fluid size and velocity, by 
the shear viscosity. High values of this ratio (low viscosities) imply turbulent, unstable 
flows, whereas small values (large viscosities) allow laminar, stationary flows. These 
ratios however are contingent on the particular flow characteristics, and not on the 
thermodynamic state of a substance alone. 

Recently, the ratio of viscosity to entropy density r\ /s in natural units, or &grj /(%s) in 
arbitrary units, has received much attention. This particular dimensionless ratio controls 
how fast shear perturbations damp, since the appropriate dissipative transverse disper- 
sion relation is co(k) = —i^jT ^ or van i snm g chemical potential). 

The ratio r\/s has been much investigated because of the work of Kovtun, Son and 
Starinets (KSS) [2], who observed that in field theories accepting a gravity dual in higher 
dimension through the holographic principle, the ratio could be estimated as a function 
of the metric coefficients near a "black brane" 



7 = Tf[g a p] 
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FIGURE 1. T7 /s for molecular nitrogen and Helium gases together with the KSS bound. 



and to their surprise, several feasible calculations with simple metrics g a p consistently 
yielded the value T]/s = 1/ (47r). The field theories dual to these gravity configurations 
are strongly coupled supersymmetric Yang-Mills theories, far removed from our current 
physical picture of the world. However, going to common substances whose viscosity 
and entropy density values are tabulated, one finds (see Fig. 1) that this ratio is at least 
an order of magnitude larger than 1 / (4n) . The reader can find details in [3] about 
simple theoretical constructions that violate this bound, but the fact remains that no 
experimental system has been so observed. 

An intuitive argument [4] for the existence of such a bound is given by Heisenberg's 
uncertainty principle, £T > 1, with x the characteristic damping time (controled as above 
by the viscosity), and £ the average energy per particle in the plasma. Some dimensional 
analysis transforms this relation in the KSS bound, excepting the 1 / (4n) prefactor. 

The ratio can also be approximately extracted in Relativistic Heavy Ion Collisions, 
and is there shown to be very close to the lower bound, which has lead to the expression 
"perfect liquid" to refer to the strongly interacting nuclear plasma (perhaps describable 
in terms of quarks and gluons instead of hadrons). The key observable is the so called 
"elliptic flow". To understand it, the reader just needs to return to Eq. (1) and start by 
looking at the ideal-fluid part of the equation. A larger gradient of pressure translates 
into a larger acceleration, while the viscosity terms damp this additional acceleration. 

In Heavy Ion Collisions the nuclei do not always collide head-on, but only part of the 
two spheres overlap. Then the pressure gradient in the reaction plane is larger than the 
pressure gradient in the perpendicular direction. The resulting momentum distribution 
of emitted hadrons peaks therefore at higher momenta for in-plane hadrons causing 
an azimuthal anisotropy around the beam axis. Its clear observation indicates that the 
viscosity is low. 



MICROSCOPIC COMPUTATIONS OF VISCOSITIES 



The one-particle out-of-equilibrium distribution function f(x,\,t) satisfies a transport 
equation 

df(x,y,t) 



dt 



or more explicitly, for a boson gas, 
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Slightly out of equilibrium, / ~ /o + 5/, where /o is the Bose-Einstein distribution 
function. The perturbation function 8f is proportional to the gradients of the hydrody- 
namical fields out of equilibrium, the velocity if one deals with viscosities. To obtain the 
shear viscosity we parametrize 8f as [5] 



f = fo 
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where Ajj = pipj — \ 5jjp 1 and Vjj represents the gradient of the velocity, v ; - : 



(8) 



Eq. (6) is a slightly involved integral equation that is solved projecting it into an 
appropriate basis of functions, performing the multidimensional integral by Montecarlo, 
and then inverting the matrix truncation of the collision operator C. 

Once solved, the shear viscosity, a macroscopic hydrodynamic coefficient, is ex- 
pressed as an integral over the distribution function, dependent on such microscopic 
details as the mass of the particles in the gas and the cross section between them, and a 
function of the temperature. 

An analogous derivation yields the bulk viscosity, except Eq. (7) is substituted by 
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and we can compare them 
1 = 



dp 

ior J E{p) 
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(having introduced the speed of sound vj. 



The computation of the bulk viscosity is complicated by the presence of zero modes 
in the corresponding collision operator C. If one substitutes gr = 1, the symmetry 
properties of C yield zero. This zero mode is associated to particle number conservation 
and is not really a departure from equilibrium. Associated to energy conservation there 
is a second zero mode E. 

Therefore one has to exclude the constant and linear oc E functions from 8f for bulk 
distortions, and invert the collision operator in the perpendicular space to these zero 
modes. 

The complications with the bulk viscosity are further enhanced by the Landau- 
Lifschitz condition, that becomes a restriction on the allowed gr, given by 

I^foE(p) 2 g^p)=0. 

SHEAR VISCOSITY AND PHASE TRANSITIONS 

It was stressed in [6] that the shear viscosity has a minimum near phase transitions in 
conventional fluids as in Fig. 1, and that this could be used in experimental work to 
pinpoint the precise location of the possible phase transition between the hadron gas and 
the quark-gluon "perfect liquid". Detailed work on the pion gas is consistent with this 
hypothesis [7, 8]. 

Since there is no formal proof of this statement, but just indications from particular 
calculations in QCD, we have studied the issue in a toy model that allows for a much 
more reliable treatment (still, in mean-field approximation for the vacuum condensate), 
the linear a model in large-/V [9]. This model is the counterexample that, although the 
phase transition and the minimum of viscosity over entropy are related, the temperature 
of both do not need to coincide. 

The situation is clearly seen in Fig. 2. The KSS number is seen to have a non- 
analiticity at the second-order phase transition where the condensate vanishes, but the 
minimum is seen to be reached distinctly before this phase transition, when the absolute 
value is dropping and affecting the cross section for pion scattering (and hence the shear 
viscosity) most greatly. 

In conclusion, although true that viscosity over entropy density is an indication of 
the phase transition, we do not dare conclude that it will pinpoint the phase transition 
temperature precisely. 

BULK VISCOSITY AND SCALE INVARIANCE 

The bulk viscosity is sensitive to the conformal breaking of the system. If the system is 
dilatationally invariant, then £ is equal to zero. If conformal symmetry is broken (e.g. 
by the explicit presence of dimensionful parameters in the Lagrangian like a mass term) 
then £ turns out to be nonzero and proportional to this breaking term. The conformal 
invariance can also be broken by pure quantum effects (trace anomaly), where the 
interaction measure (defined as the trace of the momentum-stress tensor normalized by 
the temperature T a a /T A ) is nonzero. This growth of the trace anomaly near T c has been 
proposed to drive a maximum in the bulk viscosity. 
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FIGURE 2. The minimum of rj /s occurs just before the critical temperature for the phase transition in 
the Linear Sigma Model. This is where the condensate varies rapidly approaching zero. The second order 
phase transition is marked by a non-analyticity of the ratio. 



In the pion gas the relation between the trace anomaly and the bulk viscosity has 
been presented in [10] showing a peak in £ when the effects of the trace anomaly at 
^(r 8 ) in the chiral Lagrangian are included. However, in the recent work [11] where 
the bulk viscosity in the massive Gross-Neveu model in 1 + 1 dimensions is calculated, it 
is shown that there is no clear relation between the maximum in the interaction measure 
and a peak in the bulk viscosity, showing the latter a monotonous behaviour with the 
temperature. Moreover, some of the previous results on the bulk viscosity are not always 
compatible as we can see in the Fig. 3 and show important quantitative discrepances. 

To clarify this situation we are calculating the bulk viscosity in the pion gas with 
unitarized chiral perturbation theory in the interactions both the elastic collision process 
and the speed of sound. In the low temperature limit (where the inelastic processes that 
enter are exponentially suppressed), the bulk viscosity asymptotically behaves as 



and to zero at high temperature. In the independent particle approximation [12], where 
the thermodynamics is that of an ideal gas, this occurs monotonously. The speed of 
sound has a universal behaviour in this approximation, monotonously increasing towards 
the aymptotic value of 1 / a/3. 

Our preliminary result within the independent particle approximation to Eq. 6 is 
shown in Fig. 3. The order of magnitude of £ is compatible with the previous result 
of [10], but the qualitative behaviour is strictly monotonous in agreement with [13], 
presumably because of the universality of the speed of sound. The peak around T = 230 
MeV in [10] could appear in our calculation by improving the dependence of c s on 
temperature. 



Bulk viscosity of the pion gas 
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FIGURE 3. The value of the bulk viscosity in the pion gas, as well as its temperature dependence, 
has not yet been settled. Solid line: chiral perturbation theory based calculation [10]. Dashed line: 
kinetic theory calculation based on phase shifts [13]. Dashed-dotted line: Our preliminary result in the 
independent particle approximation for the pion gas. 



We are currently investigating these discrepancies in detail, as well as lifting the in- 
dependent particle approximation including the effect of a more realistic behaviour of 
the speed of sound through the unitarized cross section coming from the Inverse Am- 
plitude Method applied to the chiral perturbation theory. The results will be published 
elsewhere. 
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